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Weather radars, recording information about precipitation around the globe, will soon be 
significantly upgraded. Most of today’s weather radars transmit and receive microwave 
energy with horizontal orientation only, but upgraded systems have the capability to send 
and receive both horizontally and vertically oriented waves. These enhanced “dual- 
polarimetric” (DP) radars peer into precipitation and provide information on the size, 
shape, phase (liquid / frozen), and concentration of the falling particles (termed 
hydrometeors). This information is valuable for improved rain rate estimates, and for 
providing data on the release and absorption of heat in the atmosphere from condensation 
and evaporation (phase changes). The heating profiles in the atmosphere influence global 
circulation, and are a vital component in studies of Earth’s changing climate. However, 
to provide the most accurate interpretation of radar data, the radar must be properly 
calibrated and data must be quality controlled (cleaned) to remove non-precipitation 
artifacts; both of which are challenging tasks for today’s weather radar. The DP 
capability maximizes performance of these procedures using properties of the observed 
precipitation. 

In a notable paper published in 2005, scientists from the Cooperative Institute for 
Mesoscale Meteorological Studies (CIMMS) at the University of Oklahoma developed a 
method to calibrate radars using statistically averaged DP measurements within light rain. 
An additional publication by one of the same scientists at the National Severe Storms 
Laboratory (NSSL) in Norman, Oklahoma introduced several techniques to perform 
quality control of radar data using DP measurements. Following their lead, the Topical 
Rainfall Measuring Mission (TRMM) Satellite Validation Office at NASA’s Goddard 
Space Flight Center has fine-tuned these methods for specific application to the weather 
radar at Kwajalein Island in the Republic of the Marshall Islands, approximately 2100 
miles southwest of Hawaii and 1400 miles east of Guam in the tropical North Pacific 
Ocean. This tropical oceanic location is important because the majority of rain, and 
therefore the majority of atmospheric heating, occurs in the tropics where limited ground- 
based radar data are available. 

The current study supports the conclusions of CIMMS and NSSL that DP radar 
measurements in light rain provide calibration results very similar to independent, but 
logistically difficult, techniques. Quantified data reveals that DP-based quality control 
results in significantly improved data quality and more accurate rainfall measurements. 
In fact, data contamination from sources such as buildings, towers, ocean waves, and 
even density changes within the air itself, can be almost completely corrected. The 
practical application of these DP methods to tropical oceanic radar data provides a 
calibrated and clean radar dataset crucial to researchers studying our changing planet. 
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Abstract 


The dual-polarization weather radar on Kwajalein Atoll in the Republic of the Marshall 
Islands (KPOL) is one of the only full-time (24/7) operational S-band dual-polarimetric (DP) 
radars in the tropics. Through the use of KPOL DP fields and disdrometer data from Kwajalein, 
quality control (QC) and reflectivity calibration techniques were developed and adapted for use 
in a near real-time operational environment. Data studies in light rain show that KPOL DP 
measurements meet or exceed quality thresholds for these applications as determined by 
consensus of the radar community. While the methodology for development of such applications 
is well documented, tuning of specific algorithms to the particular regime and observed raindrop 
size distributions requires a comprehensive testing and adjustment period. Presented are 
algorithm descriptions and results from five case studies in which QC and absolute reflectivity 
calibration were performed. Also described is a unique approach to calibrate the differential 
reflectivity field when vertically pointing scans arc not available. Results show the following: 1) 
DP-based QC provides superior results compared to the legacy Tropical Rainfall Measuring 
Mission (TRMM) QC algorithm (based on height and reflectivity thresholds), 2) absolute 
reflectivity calibration can be performed using observations of light rain via a published 
differential phase-based integration technique; results are within ± 1 dB compared to 
independent measurements, and 3) a polarimetrically tuned reflectivity-rain rate (Z-R) 
application extended to six months of data shows very good agreement with independent rain 
gauge measurements, thereby extolling the benefit of DP-based QC and calibration of the 
reflectivity field. 
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1. Introduction and KPOL site description 

Dual-polarimetric (DP) ground-based weather radars are well recognized as vital 
instruments for applications in hydrology, precipitation microphysics, and hydrometeor 
identification (Ryzhkov and Zmic 1998a, Vivekanandan et al. 1999, Straka et al. 2000, Gorgucci 
et al. 2001, Wang and Carey. 2005, among others); all of which benefit research focusing on 
ground validation of satellite measurements (Chandrasekar et al. 2008). Kwajalein, Republic of 
the Marshall Islands (KWAJ) (Figure 1) is an ideal tropical oceanic location for which ground 
validation, modeling, mesoscale characterization studies, and other activities have focused 
(Schumacher and Houze 2000, Sobel et al. 2004, Houze et al. 2004, Wolff et al. 2005, Yuter et 
al. 2005, Blossey et al. 2007, Wang and Wolff, 2009). In support of U.S. Army and Tropical 
Rainfall Measuring Mission (TRMM) Ground Validation (GV) operations at KWAJ, an S-band 
DP radar (KPOL) operates on a continual basis, providing unique opportunities for operational 
algorithm development and adaptation with applications clearly extendable to the Global 
Precipitation Measurement (GPM) GV program. A goal of the GPM GV program is to improve 
accuracy of rainfall retrievals by developing and improving physically based radiometer 
algorithms for application over land and ocean. This approach requires insight into the 
properties of ice microphysics, parameters of local and regional drop size distributions (DSDs), 
and delineation between water phases. Through careful observation, DP radars can provide a 
means to cross-validate parameterized microphysical properties in GPM radiometer retrieval 
algorithms (Kummerow and Petersen 2006, Chandrasekar et al. 2008). As explained in these 
references, the utility of any operating radar and rain gauge networks for establishing random 
and systematic errors in rainfall relative to the GPM core satellite over oceanic domains should 
be fully explored. Kwajalein is an oceanic site, and the operational DP algorithms for quality 
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control (QC), self-consistency calibration, hydrometeor identification, and rainrate estimation 
provide an opportunity for validation of microphysical properties in ocean-based radiometer 
retrievals from storm to climate scale. Over land, cloud resolving models (CRMs) are one of 
several components that will be utilized in the development of ground-based physical retrieval 
algorithms (Kummerow and Petersen 2006). The operational DP algorithms can be extended to 
WSR-88D radar sites thereby providing needed input to CRMs with currently parameterized 
microphysical properties. Deployment of DP diverse WSR-88D radars is scheduled to begin in 
October 2010 (Istok et al. 2009). 

The primary focus of this manuscript is the development and adaptation of operational 
algorithms for QC and self-consistency calibration with the KPOL radar at Kwajalein for the 
purpose of exploring DP-based validation capabilities. Section 2 describes the overall quality of 
the KPOL data as compared to established DP research radars, and details the physically based 
QC techniques. Section 3 provides a method for calibration of differential reflectivity, and 
application of self-consistency reflectivity calibration using properties of the rainfall medium. 
DP calibration results are compared with an independent statistical calibration approach. Finally, 
the benefit of a DP-based QC’d and calibrated reflectivity field is demonstrated with application 
of the polarimetrically-tuned Z-R method of Bringi et al. (2004). 

2. KPOL Quality Control 

a. Data quality 

KPOL was upgraded to DP capability in early 1998. Routinely observed parameters and 
scanning strategies are shown in Tables 1 and 2, respectively. Due to a litany of engineering and 
mechanical issues with KPOL (Marks et al. 2009), initial DP data were not up to established 
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standards of acceptability until 2006. To determine the basic quality of KPOL DP data, 
empirical comparisons were made with established DP research radars. The National Center for 
Atmospheric Research (NCAR) S-band polarimetric radar (SPOL) (Lutz et al. 1995) and the 
Colorado State University - University of Chicago - Illinois State Water Survey (CSU-CHILL) 
(Brunkow et al. 2000) radar were used as comparison benchmarks to evaluate relative KPOL 
performance. This was accomplished through analysis of DP measurements in very light rain 
(20 dBZ < Zh < 28 dBZ, where Zh is the horizontal reflectivity component). In this context, 
drops are essentially spherical with little or no variability in shape, canting angle, or scattering 
properties within a radar resolution volume (Doviak and Zmic 1993). Therefore, DP 
measurements from this medium are used as indicators of data quality. In polarimetric radar 
applications, the co-polar correlation coefficient (phv) is a measure of the correlation between 
horizontally and vertically polarized weather signals. It is primarily affected by the variability in 
the ratio of the vertical-to-horizontal size of hydrometeors in the resolution volume, but can also 
be affected by variability in canting angle and differential phase shift on scattering (Doviak and 
Zrnic 1993). Canting and scattering properties are of minimal concern in light rain and drizzle 
due to the mainly spherical shape of the small drops. Theoretical values of pnv larger than 0.99 
are indicated by Sachidananda and Zrnic (1985) because of the small shape effects in rain, but 
theory does not account for possible decreases due to sidelobes and receiver noise. 
Measurements in rain indicate an average pnv of 0.98 (Doviak and Zmic 1993) with standard 
deviation of 0.01. Therefore, significant deviation (> 0.01) of pnv below 0.98 in light rain is a 
likely indicator of general radar system issues. Other values of polarimetric measurements in 
light rain (as shown in Doviak and Zrnic 1993 Tabic 8.1) include the analysis of both specific 
differential phase (Kdp) and differential reflectivity (Zdr), where median Kdp measurements 
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should be approximately 0° km" 1 , and deviation of Zdr from its average value should be 
approximately 0 dB, i.e. 

£(Z Dfi (0-Z^)/A^0dB (1) 

l'=l 

where N is the Zdr sample size. DP observations in light rain were compared from the three S- 
band radars, SPOL, CHILL, and KPOL. Figure 2 shows representative pnv, deviation of Zdr, 
and K D p measurements by relative frequency for specific events from the three radars under 
typical operating conditions. Base tilt data from representative volume scans were compared. 
All three panels of Figure 2 indicate that KPOL data quality is at least as good as the established 
research radars in this specific case. To investigate further, twelve KPOL volume scans from 
different light rain events from 2006 and 2007 were analyzed with very similar results. To 
determine if Kdp data are of sufficient quality for applications of QC and rain rate estimation, 
phase data were also examined. The standard deviation of measured differential phase, o(Odp), 
was computed at each range gate from a running, centered 25-gate sample in the radial direction. 
This corresponds with 0.2 km gate spacing and the RVP8 Kdp length scale of 5 km. As 
discussed in the literature, a reasonable range for g(Odp) computations should be approximately 
2-3° or lower for Kdp applications involving QC, calibration, rain rate estimation and 
hydrometeor identification (e.g., see Doviak and Zmic 1993, and Bringi and Chandrasekar 2001 
for review). In addition, the absolute average deviation (AAD) of Zdr from its mean should 
provide an upper-end estimate of the random measurement error in Zdr. The AAD of Zdr is 
defined by equation 2, with N being the number of Z D r measurements. 

AAD[Z dr ]=X|(Z m (/)-Z^)|W (2) 

1=1 
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Table 3 shows statistical relationships between median Kdp, pnv, cj(Odp), and AAD of Zdr from 
the CHILL, SPOL, and KPOL radars in light rain. The KPOL statistics represent the range of 
values from twelve analyzed light rain cases from 2006-2007. These results show that KPOL 
data quality meets and typically exceeds the radar community’s consensus data standards for the 
implementation of DP radar-based QC, self-consistency calibration, rain rate estimation, and 
hydrometeor identification algorithms. 
b. Radar Data and QC techniques 

KPOL reflectivity and DP data are frequently contaminated by ground and sea clutter, 
multiple-trip echo, and considerable noise. QC algorithms based on DP measurements have 
shown notable success in objective identification of these and other non-precipitation features 
(Ryzhkov and Zmic 1998b, Zmic and Ryzhkov 1999, Cifelli et al. 2002). A series of five KPOL 
case studies with varying rainfall intensities and reflectivity coverage from years 2006 and 2007 
were selected to develop an operational QC algorithm for detection and removal of non- 
precipitation echo. A monthly study from July 2008 was also performed to check the DP QC 
algorithm in various other conditions such as isolated to scattered convection and light showers, 
to periods with only non-precipitation echo present. Since the primary focus was in detection 
and removal of non-precipitation echo during rain events, the discussion concentrates on results 
from the five individual case studies. In the following discussion, refer to Table 1 for description 
of moments and field labels. The multiple steps of the QC algorithm for the first two elevation 
angles are displayed in flowchart form in Figure 3. For a primary application of quantitative 
rainfall estimation, the QC steps are applied below the observed radar bright band and site- 
specific melting level of approximately 5-6 km (Schumacher and Houze, 2000). In our DP QC 
algorithm, a new data field with label “CZ” is created for each volume scan and contains the 
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final pre-calibrated reflectivity that has been edited for suspected non-precipitation echo. 
Initially, the CZ field is simply copied from the raw reflectivity field (ZT) for all elevations scans 
within a volume. As gates are identified as non-precipitation echo, they are assigned a specific 
value corresponding to the no-data flag. The first step in the QC process is automated and 
applied by the RVP8 processor. A signal-to-noise ratio (SNR) test is applied to the Odp(PH), 
Kdp(KD), Zdr(DR), and pnv(RH) fields to identify gates with weak or uncertain signals, and sets 
the value of these gates to the no-data flag. Multiple-trip echo is usually removed from the Odp 
and pHv fields by this technique; however, additional QC is required for multiple-trip in both Zh 
and Zdr fields. All gates containing the no-data flag in the Odp, Phv, Kdp, and Zdr fields are 
mapped to the corresponding gates in the QC’d reflectivity field CZ. This step takes advantage 
of the relatively clean Odp and pnv fields, and eliminates those gates from CZ that have been 
flagged with low SNR by the RVP8 processor. This step also removes gates from CZ for which 
there are no corresponding DP measurements. 

The calculation of o(Odp) at each range location, and subsequent threshold comparison, is 
shown in Ryzhkov and Zmic (1998b) to be a successful test for detection of anomalous 
propagation (AP)-induced ground clutter echo. The technique is also effective for any type of 
ground clutter or region of low SNR, and has been applied to KPOL data for identification of 
ground clutter associated with human-made structures (i.e. buildings and towers). Before 
computing the standard deviation, the Odp field is de-aliased. The total differential phase 
contains both the radar system phase and a cumulative phase due to scattering from precipitation 
(Gorgucci et al. 1999; Bringi and Chandrasekar 2001). KPOL system phase has a history of 
variation from near 0° to just shy of the maximum unambiguous value of 180°. As a 
consequence, the Odp field has shown varying amounts of aliasing depending upon the data 
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being analyzed. A gate is considered aliased if the absolute value of the phase difference 
between consecutive gates exceeds 149°. An additional correlation test is performed to make 
certain that noise is not being sampled. All aliased gates are corrected by adding 180° to the 
existing phase value. Using a running centered 15-gate sample, o(Odp) is then computed at each 
range gate. If at least 5 of the 15 gates contain valid Odp measurements, their standard deviation 
value is assigned to the center of the radial interval; otherwise the o value is set to the no-data 
flag. The requirement of five or more phase samples for a standard deviation calculation 
eliminates isolated speckle. When o(®dp) calculations are complete for a given sweep, each gate 
value is checked against an empirically determined threshold. Within precipitation (15 dBZ < 
Zh<55 dBZ) the typically observed average o(Odp) from KPOL is about 3°. Choosing a multiple 
of four, the threshold has been set to 12° and is the same threshold used in Ryzhkov and Zmic 
(1998b). 

Sea clutter, ground clutter from structures (buildings and towers), and general noise are 
detected and eliminated by using a combination of o(O dp ) and pnv thresholding. If a o(O dp ) gate 
is greater than the threshold, or has been set to the no-data flag, the corresponding CZ gate is set 
to the no-data flag. Similarly with pnv, if a correlation gate is less than the threshold of 0.80, the 
corresponding CZ gate is set to the no-data flag. Analysis of p H v within sea clutter reveals 
values mostly less than 0.40; however values in the range from 0.0 to near 0.95 can occur. A 
similar analysis of o(Odp) within sea clutter shows standard deviation values ranging from 3° to 
70°. This wide range of correlation and standard deviation values is expected given the varying 
nature of returns from ocean waves. With a pnv threshold of 0.80, and a g(Odp) threshold of 12°, 
almost all sea clutter is detected. Echo clearly identified as ground clutter from reflectivity time- 
series analysis displays typical pnv values in the 0.4 to 0.95 range. More than 50% of these 
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ground targets have pnv values exceeding 0.80, and could easily be incorrectly identified as 
precipitation echo if the correlation test was considered alone, therefore, the o(0 dp) test is also 
needed. Within ground clutter, o(Odp) has values ranging from 10° to near 80°, with a clear 
majority of values greater than 40°. The combination of the correlation and standard deviation 
tests identifies almost all ground clutter gates; however a small percentage of problem gates are 
not flagged by either threshold and survive the QC tests. 

An additional series of QC steps to test for realistic values and noise are applied to the 
Kdp and Zdr fields. The Zdr field is calibration corrected prior to this test, and is the topic of the 
next section. The specific thresholds applied in these tests are appropriate for KPOL, but can 
easily be tailored to match the nature of site-specific data at other locations. Even in the heaviest 
rain at Kwajalein (Zh near 52 dBZ), Kdp values do not exceed 2.4°km" 1 as determined from 
analysis of both radar and disdrometer data. The K D p realistic value test considers gates with 
values > 3°km’' or < -2°km"' as bad data or noise. The positive end of the Kdp allowed range 
corresponds to a rain rate of -136 mm hr" 1 using the R(Kdp) equation of Bringi and Chandrasekar 
(2001). Analyses of Zdr by reflectivity bin (post calibration) from the lowest two elevation 
angles (0.4° and 1.4°) indicate average values approaching 1-1.25 dB for the heaviest rain, and is 
a similar Z D r range computed from disdrometer observations (discussed in the following 
section). For a primary application of quantitative rainfall estimation from the lowest angles, the 
Zdr noise test considers gates with negative values, or those > 2.5 dB to be suspect. The upper 
threshold was subjectively determined based on low-elevation KPOL observations under typical 
operating conditions. These additional tests are needed as spurious Kdp and Zdr gates can be 
persistent in contaminating the interpolated fields for subsequent rain rate estimation. The 
percentage of gates flagged as “bad” by each of the QC tests from the first two elevation angles 
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and rainfall estimation range of 0-100 km is shown in Table 4. Within these elevation and range 
limits, we are well below the melting level and radar observed bright band. The higher 
percentage of flagged gates for the Zdr threshold test is due to values less than 0 dB. Gates 
flagged in total reflectivity (ZT) includes those from the pnv, <j(Odp), low reflectivity threshold 
(< 5 dBZ), and missing gate tests from the DP fields. The ZT percentages vary significantly 
from case-to-case and are heavily influenced by flagging those gates less than 5 dBZ. 

The CZ field now contains the corrected (pre-calibrated) data. As a final QC step for 
noise reduction in DP fields, all CZ gates containing the no-data flag arc mapped to the 
corresponding gates in pnv, Zdr, Kdp, and Odp. The first two elevation sweeps are QC’d in this 
manner. Figure 4 shows typical results of the QC algorithm for the 0.4° Zh field. The top panels 
(a,b) of Figure 4 show raw and corrected reflectivity images within 50 km radius and indicate the 
effective identification of ground clutter along the atoll perimeter (both embedded and non- 
embedded) in precipitation echo. It is clear that sea clutter, multiple-trip echo, and general noise 
have also been identified and removed. The bottom panels (c,d) of Figure 4 show a full 160 km 
0.4° sweep before and after QC. Pronounced regions of multiple-trip echo (from 220 deg to 250 
deg) have been removed in addition to widespread light noise. 

A quantitative analysis of ground clutter returns reveals that DP QC results are superior 
to the legacy TRMM GVS (Ground Validation System non-polarimetric) algorithm. As 
discussed in Kulie et al. 1999, the GVS QC algorithm identifies non-precipitation echo by use of 
height and reflectivity threshold parameters, and has a significant weakness in removal of high- 
reflectivity ground clutter, especially when the clutter is near or embedded within precipitation. 
The strongest precipitation echoes at Kwajalein approach 52 dBZ, but ground clutter returns 
easily exceed this value with measurements ranging from 55 to 70 dBZ. Figure 5 shows the 


11 



location of the clutter field at Kwajalein, with 1323 gates (within 50 km of KPOL) identified as 
frequent sources of clutter (Silberstein et al. 2008). Reflectivity gates are extracted exclusively 
from these locations from unedited (raw), and corrected data from both DP and GVS QC 
algorithms for the five daily case studies. To be reasonably certain that no precipitation echo is 
selected, only reflectivity values > 55 dBZ are considered to be ground clutter. Table 5 shows 
the gate count > 55 dBZ from uncorrected, DP corrected, and GVS corrected data. The 19 Dec 
2006 case shows that 71 gates from a total of 1 1907 extracted gates have values > 55 dBZ. DP 
QC has correctly identified and removed all 71 gates (100% correction), therefore zero clutter 
gates remain. GVS QC has 53 remaining clutter gates > 55 dBZ, roughly corresponding to a 
25% correction. Similar results are shown for all cases. DP QC has virtually no clutter gates 
remaining, while GVS QC has significant numbers of clutter gates remaining. In these cases, 
precipitation echo is widespread and covers significant portions of the field, and ground clutter 
echo is mostly embedded in (or in close proximity to) precipitation echo. The DP QC tests 
(correlation and standard deviation of phase) detect and remove the embedded clutter, but GVS 
QC historically fails in this regard. In cases with partial precipitation coverage and non- 
embedded clutter, it is possible for marginal improvement of GVS QC performance through 
threshold strengthening, but requires repetitive labor intensive processing. In contrast, the DP 
QC algorithm is fully automated and provides consistent results without the requirement of 
parameter adjustments. 

3. KPOL calibration 

a. Differential reflectivity calibration 

Accurate Zdr calibration is essential in the determination of absolute reflectivity 
calibration (Zbias) via consistency among the polarimetric variables Zh, Zdr, and Kdp, and is 
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also critical for rain rate estimation and hydromctcor identification. Approaches to calibrate the 
differential reflectivity field can be both engineering oriented and through the use of natural 
scatterers. The use of vertically pointing (or birdbath) scans in light rain is a favored and reliable 
approach to determine Zdr bias (Gorgucci et al 1999, Hubbert et al 2008). The KPOL dataset 
from years 2006 and 2007 does not contain reliable birdbath scans; therefore an alternative 
calibration method was needed. Zmic et al. (2006) describe an approach to partition Zdr 
calibration into measurements of constant system biases and time-varying biases with the 
assumption that the time constant of any slow calibration drift is substantially longer than the 
duration of a volume scan. Their method is applicable to WSR-88D systems due to the need for 
calibration results without reliance on vertical profiles in precipitation. This approach is not 
easily applicable to KPOL without hardware modifications. Hubbert et al. (2003) perform Zdr 
calibration by employing a cross-polar power technique from solar measurements. In their 
approach, solar scans are used to calibrate the receiving path, and the average cross-polar powers 
for horizontal and vertical orientations are obtained from the backscatter covariance matrix to 
calibrate the transmission path. However, this technique is not applicable to KPOL due to the 
unavailability of cross-polar power measurements while operating in simultaneous transmit and 
receive mode. Using natural properties of scatterers, Ryzhkov et al (2005) analyzed dry 
aggregates at high elevation angles above the radar-detected melting level. The dry aggregate 
medium has very low density and therefore has small intrinsic Zdr measurements. This method 
may be applicable for individual case studies, however, for the operational environment of 
KPOL the dry aggregate approach was found to yield inconsistent results from scan-to-scan. In 
practice, it is very difficult to isolate the data above the melting level, avoid contamination from 
the melting level itself, and yet not be sampling ice medium with preferred orientations. Bcchini 
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ct al. (2008) propose a calibration method based on properties of the rain medium where Zdr 
measurements are collected at increasing elevations. The Zdr profiles are then compared to 
theoretical profiles to detect and quantify bias. Their method was applied to an operational radar 
where the scanning strategy prohibited vertical observations. 

Instead of developing a theoretical Zdr profile as in Bechini et al. (2008), a unique 
method was applied where the calibration of KPOL Zdr was accomplished through bias 
adjustment to a disdrometer-detennined Zdr reference profile. Over 10 000 JW disdrometer 
observations of Zh and Zdr from 2003 and 2004 were compiled for the reference profile. 
Assumptions regarding drop size and shape relations used in the disdrometer Zdr computation 
are described in the following section. KPOL Zdr data from 2006 and 2007 smoothed via a 
boxcar approach (Bringi and Thurai, 2008) were then calibrated for individual rain events by 
application of specific offsets as determined by comparison to the disdrometer reference. The 
Zdr Zh disdrometer reference profile is shown in the top panel of Figure 6 (bold line with no 
symbols) together with profiles from five case studies in 2006 and 2007. The cases were chosen 
based upon rainfall coverage, and include those with uniform rain shields containing small 
embedded convective cells. The level of disagreement in Zdr distributions within the cases is 
evident, and their bias relative to the disdrometer reference is shown. Before determining the 
proper Zdr offset, the Zh distributions were independently calibrated by the relative calibration 
adjustment (RCA) method (Silberstein et al. 2008; Marks et al. 2009). The bottom panel of 
Figure 6 shows the Zdr distributions after adjustment to the reference. Emphasis was placed on 
matching within the 30-40 dBZ range as this represents approximately 85-90% of the 
measurements used in self-consistency calibration (section 3b). A trial-and-error adjustment of 
Zdr offsets with comparison to the disdrometer reference resulted in the final offsets to apply in 
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each case. While the adjusted Zdr profiles are not in perfect agreement, and fluctuations due to 
sample size limitations are noticeable especially in the upper reflectivity bins, our analysis 
indicates that the agreement is sufficient to perform a robust calibration of Zh using the self- 
consistency technique. 

Occasional performance of birdbath scans in light rain commenced in March 2008, and 
have continued through 2009. Results consistently show Zdr running ‘hot’ by approximately 1.0 
dB. Figure 7 displays birdbath scan data from 07 Mar 2008 during a shallow light rain event. 
The vertical profile of mean Zdr (top right panel) shows bias of approximately +1 dB at the 3.0 
km level; a height chosen to be above pronounced ground clutter effects and below the melting 
level. The mean pnv values near 0.99 are indicative of sampling in the rain medium (lower left 
panel). A 3.0 km mean Zdr distribution by azimuth in the bottom right panel shows the periodic 
nature of the direction-dependent measurements; an expected structure due to variability in 
ground clutter response with antenna rotation in the vertical (Gorgucci et al. 1999). Zdr is 
averaged from the full azimuth cycle of 360°, over which 1023 samples per azimuth are 
obtained. 

b. Self-consistency calibration: methodology ? and results 
i. Methodology 

Absolute radar calibration is a requirement for quantitative rainfall estimation (Ulbrich 
and Lee 1999). A calibration offset of 2-dB (for example, from 28 to 30 dBZ) can result in 
rainfall estimation error of -30% as measured by the default WSR-88D reflectivity rain rate 
relationship (Houze et al. 2004). Some early and relatively simple calibration approaches are 
described in Atlas (2002). It is well documented that polarimetric properties of the rain medium 
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can be used to determine the absolute calibration of a radar system. Techniques to capitalize on 
these consistency relations range from the comparison of rainfall rates derived from power and 
phase measurements (Gorgucci et al. 1992), to comparing observed and estimated differential 
propagation phase (Goddard et al. 1994, Vivekanandan et al. 2003, Ryzhkov et al. 2005, and 
others). The self-consistency of Zh, Zdr, and Kdp measurements was quantified by Scarchilli et 
al. (1996) and Gorgucci et al. (1999) using a gamma distribution model that described many of 
the natural variations in the DSD. The majority of these self-consistency calibration techniques 
have a common necessity of heavy rain rates (> 50 mm h" 1 ) for significant phase accumulation at 
S-band over the range profile. Ryzhkov et al. (2005) developed consistency relationships based 
on existing statistics of DSD measurements and polarimetric radar observations in central 
Oklahoma, and suggested a methodology for determining absolute reflectivity bias (Zbias) from 
the self-consistency relation that did not require heavy rainfall at S-band. This methodology 
compared area-time integrals of measured (processor or user determined) Kdp and computed 
(theoretical) Kdp (as a function of Zh and Zdr) and determined Zbias as the adjustment in Zh 
needed for the integrals to agree. 


Integrating Kdp over a large space-time domain substantially reduces the inherent 
noisiness of point K D p measurements, thereby allowing lighter rains with relatively low <D D p to 
be acceptable for self-consistency calibration of Zh. Per Ryzhkov et al. (2005), to obtain a self- 
consistency Zbias within 1-dB, the area-time integral of measured Kdp should be estimated 
within 20% accuracy. To determine the accuracy or standard deviation of Kdp, we use the 
expression from Balakrishnan and Zmic (1989) and Carey et al. (2000) (equation Al) 


ct(K dp ) = 


\[3ct(<S> dp ) 

N L5 A r 


( 3 ), 
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where a(0 DP ) is the standard deviation of the total differential phase, N is the number of range 
gates used in Kdp calculations, and A r is the range gate spacing. The length scale for KPOL Kdp 
calculations is 5-km and the gate spacing is 0.2 km, therefore, a 25 gate filter is used for each 
Kdp calculation. For the five case studies, o(Odp) was calculated over 25 gate blocks from all 
volume scans with careful attention to include data within specified range and azimuth 
boundaries that contained continuous echo (i.e. phase data from echo edges can be excessively 
noisy, therefore, these observations are avoided in the standard deviation calculations). The 
lowest o(Odp) observed was 2.40° from the 19 Dec. case, and the highest was 3.42° from the 1 1 
Aug. 2007 case. Inserting the lowest and highest o(Odp) values in (3), the resulting range of 
o(Kdp) in the data is from 0.17 °km‘ 1 to 0.24 °km" 1 . As part of the calibration method, Kdp 
observations are averaged within each reflectivity bin, and there is an associated reduction in 
o(K D p) relative to the number of statistically independent samples. The total number of K D p 
samples within each reflectivity bin divided by 25 (block averaging window) provides the 
number of statistically independent Kdp samples. As explained in Ryzhkov et al. (2005), the 
times reduction in o(K D p) due to averaging is defined as the square root of the number of 
independent samples. The times reduction is then applied to the o(K D p) high and low values to 
determine the o(K DP ) range in each reflectivity bin. The resulting range of a(K D p) is then 
compared to the maximum allowed o(Kdp) to assess data validity. 

This statistical technique was applied to each reflectivity bin from Z mm (30 dBZ) to Z max 
(48 dBZ) for all the case studies. The results from the 19 Dec 2006 case are shown in Table 6 
and are comparable to the results from all cases. For each reflectivity bin, Table 6 shows the 
average Kdp, the number of independent Kdp samples, the times reduction in o(Kdp) due to 
averaging, the maximum allowed o(Kdp), and the o(K D p) observed (after reduction from 
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averaging). To further clarify using the 30 dBZ bin in Table 6 as an example, the average Kdp is 
0.053°km“ 1 . There were 43526 total Kdp samples, so the number of independent samples is 1741 
(43526/25), resulting in a 41 (1741 05 ) times reduction in o(Kdp). The range for g(Kdp) becomes 
0.004 °km' 1 (0.17/41) to 0.006 °km" 1 (0.24/41). The maximum allowed o(Kdp) is 0.011 °km" 1 
(20% of the average Kdp). In this reflectivity bin, and in all bins from all cases, it is shown that 
the observed range of o(Kdp) is below the maximum allowed o(Kdp). Due to extensive 
averaging, the Kdp data are valid for use in the presented self-consistency calibration method. 

The precipitation at Kwajalein is dominated by systems that form in the intertropical 
convergence zone (ITCZ), and shallow (< 5 km) “warm rain” clouds (Schumacher and Houze 
2000, Wolff et al 2005) and is ideal for self-consistency calibration due to this lighter rain 
regime. To apply the area-integration methodology to KPOL data, DSD measurements from the 
Kwajalein region were required to derive a consistency equation between Z H , Z D r, and K D p. 
Using simulated DSD, Vivekanandan et al. (2003), for example, derived a relationship where Kd P 
is expressed as a function of Zh and Zdr. In this study, we derived a similar relation using actual 
disdrometer observations. An impact type Joss-Waldvogel (Joss and Waldvogel 1967) 
disdrometer sited at Kwajalein from May through December 2003 provided 8779 1 -minute 
resolution DSD measurements (within the 30-48 dBZ interval) to regress the following relation 
between the variables, 

K dp = AZ b H Z c DR (4) 

(where A = 0.17737xl0" 4 ; b=0.9926, and c = -0.5138) with Zh in mm 6 m°, Zdr in dB, and Kdp in 
deg km" 1 . The polarimetric radar parameters of Zh, Zdr, and Kdp, were calculated for each 
minute of DSD observations for an S-band radar (10.7 cm) and a temperature of 20 °C as shown 
in Tokay et al. (2002). For drop shape, the mean axis ratios offered by Andsager et al. (1999) 
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were adopted for drops less than 4 mm in diameter and equilibrium drop shapes (Beard and 
Chuang 1987) for larger drops. For the fall velocity, we adopted the terminal fall velocity drop 
diameter relation given by Beard (1976). The coefficient and exponents were derived via a 
linear least-squares fit regression. 

In Ryzhkov et al (2005), their derived consistency equation from DSDs in Oklahoma 
reduced the impact of variability in the DSD and raindrop shape on the calibration results due to 
the large-scale integration technique. Following their method, we matched measured Kdp and 
computed Kdp(Zh, Zdr) by adjusting Zh by an amount (in dB) considered as the Zbias- A 
practical approach to accomplish this is to divide the data collected from an entire 
spatial/temporal domain into 1-dB increments of radar reflectivity and compute average values 
of Kdp(Z) and Zdr(Z) in each 1-dB interval of Z between Z m i n (30 dBZ) and Z max (48 dBZ). The 
Zbias is them determined by matching the integrals 

^max 

/,= j K dp {Z) n(Z)dZ (5) 

'-min 

and 


max 

I 2 = A j ’z b H Z^n{Z)dz 


with an estimated Zbias determined from Vivekanandan et al. 2003 by 

Z hias (dB) = lOlogf 




( 6 ) 


( 7 ) 


An iterative adjustment approach is required to force agreement of the integrals to within an 
established bound of 0.1 dB. In all cases analyzed, this has been accomplished with two or less 
iterations. The same cases analyzed for QC and Zdr calibration were examined for self- 
consistency calibration. 


19 



ii Sclf-consistcncy calibration: Results and method verification 
As a basis for evaluation, the self-consistency results arc compared against those from the 
independent RCA technique (Silberstein et al. 2008). The RCA uses a statistical ensemble of 
reflectivity values from persistent ground clutter areas from every volume scan to monitor hourly 
and daily radar sensitivity changes relative to an established baseline. As detailed in Silberstein 
et al. (2008), the 95 lh percentile of the clutter area reflectivity distribution at the lowest elevation 
scan is remarkably stable to within ± 0.5 dB, and therefore permits monitoring of radar stability. 
A practical application of the RCA technique to KPOL data revealed a dramatic improvement in 
KPOL data stability as evidenced by reflectivity comparison with the TRMM Precipitation Radar 
(PR). Although the RCA provides a relative calibration, corrected KPOL reflectivity matched 
the PR to within ± 1 dB on a monthly and yearly basis (Marks et al. 2009; Wang and Wolff, 
2009). Table 7 shows self-consistency calibration results as compared to the RCA approach and 
the absolute value of their difference. From the case studies analyzed, there is agreement 
between self-consistency and RCA to within ±1 dB. In four cases, the agreement is within 0.5 
dB. These results arc similar to Ryzhkov et al. (2005) upon comparison of corrected reflectivity 
with independent measurements. Illingworth and Blackman (2002), and Vivekanandan et al. 
(2003) have also determined that the accuracy of self-consistency calibration can be as good as 
0.5-1 dB as evaluated through independent comparison. The KPOL results are consistent with 
previous calibration studies, and provide confidence in the operational method. As explained in 
Silberstein et al. (2008), there is ± 0.5 dB uncertainty in RCA measurements. Together with the 
possible ±1 dB uncertainty shown here, there is a combined calibration uncertainty of ± 1.5 dB. 
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To avoid spurious values of Kdp due to adverse effects from non-uniform beam filling 
(NBF) and have the best possible dataset for Kdp comparison, cross-beam gradients of Odp 
should be avoided. NBF may cause perturbations in Odp (Ryzhkov and Zrnic 1998c; Gosset 
2004), subsequently resulting in spuriously large Kdp (of both signs) in regions of strongly non- 
uniform precipitation (Ryzhkov 2007). However, for use in an operational environment, data 
within the range of 20-80 km from KPOL, and full sweep azimuths (0° - 360°) were selected 
with the belief that questionable Kdp values near echo edges would be statistically insignificant 
to the final calibration result. As suspected, due to the large amount of integrated data in each 
case study, the final calibration results using specific range/azimuth thresholding compared to 
the operational settings were very similar (within 0.25 dB). 

The consistency relationship given in (4) was derived from disdrometer observations at 
Kwajalein, and therefore the relationship can serve as a theoretical model for verification of the 
self-consistency method after bias correction. Figure 8 shows comparisons between average Kdp 
measurements (after Zbias correction) and those calculated from self-consistency [Kdp(Zh,Zdr)]. 
There is generally good agreement between the profiles in most cases (with the 23 Nov. 2007 
case being an exception). A noticeable disagreement in profiles is apparent at both lower and 
upper ends of the reflectivity range. From approximately 30-35 dBZ, SIGMET calculated K D p is 
higher than consistency theory (Zh, Zdr model). This could be due to difficulty in extracting the 
true Kdp signal from the embedded noise field despite substantial sample size and extensive 
averaging. The agreement is best within the 35-43 dBZ range, attributed to a stronger Kdp signal 
and good sample sizes. In reflectivity bins greater than 43 dBZ, the radar measured Kdp falls 
lower than the model in all cases. The number of independent Kdp samples in reflectivity bins 
from 43-48 dBZ is significantly lower than in the other bins (< 500 from the 19 Dec case - Table 
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6). Probability distribution functions show that in all cases, approximately 95% of the Kdp 
values within the distribution are accounted for within the 30-43 dBZ range. This result is not 
surprising due to the predominant occurrence of lower reflectivity precipitation at Kwajalein 
(Schumacher and Houze 2000, Wolff et al. 2005). These calibration results are obviously 
weighted with the majority of the samples; therefore not much significance can be placed on the 
results from the higher reflectivity bins. Other explanations for disagreement in profiles at the 
lower and upper ends of the reflectivity range could be related to possible bias in the 
disdrometer-derived consistency equation, and effects of processor filtering of differential phase 
for Kdp calculation. The relatively poor profile agreement in the 23 Nov. 2007 case is likely 
related to the nature of the precipitation itself. As compared to the other cases, there is 
significantly more convection present and the echo shield displays less of a uniform coverage 
pattern. The most consistent results occur when cases with uniform rain shields arc chosen. The 
calibration results presented in Table 7 and Figure 8 provide further confirmation of the Ryzhkov 
et al (2005) method, and show that an absolute bias adjustment to Zh can be determined by 
matching Kdp profiles in relatively light rain, provided a local consistency relationship is 
available. 

4. Application 

A demonstrated application of the benefit of QC’d and calibrated data is presented by 
estimating rain rate via the polarimetrically-tuned Z-R approach (PTZR) of Bringi et al. (2004), 
with six months of both corrected and uncorrected radar data. As background reference, the 
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Bringi ct al. (2004) PTZR approach begins with the assumption of a first-guess Z-R relationship 
of the form: 


Z = aR 15 (8 ) 

The a coefficient in (8) is adjusted on a pixel-by-pixel basis as the DSD and observed DP values 
evolve in space and time. A normalized gamma DSD is assumed (Testud et al. 2001), with the 
shape parameter (//), median drop diameter (Do), and concentration (N w ) obtained via S-band 
polarimetric radar measurements of Zh, Zdr, and Kdp in different rain rate conditions following 
the method of Gorgucci et al. (2002). The reflectivity factor Z is then related to rain rate R by 


z =^/? 15 




(9), 


where a'(ju) depends on // and other constants (derivations in Appendix of Bringi et al. 2004). 
Finally, they define a new a coefficient to replace the coefficient in (8), as: 


N 


(10) 


The new coefficient is calculated at each pixel, and the application of the PTZR to a 
given pixel is dynamically determined by the magnitude of observed Z H , Z DR , and K DP . This 
technique has shown very encouraging results at Kwajalein with QC’d and calibrated data. For 
example, after application of the PTZR to six months of corrected KPOL data (July-December 
2008), the resultant radar accumulations showed very good agreement with fully independent 
rain gauges, and indicate an absolute bias of only 16.5% (Figure 9 - left panel). The results from 
uncorrccted data are also shown in Figure 9 and indicate an absolute bias with rain gauges of 
74.5%. The benefit of calibrated and QC’d data is clearly evident. For this example, the 
corrected Zh data have RCA calibration adjustments applied on a daily basis, following the 
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method detailed in Marks et al. (2009). What is especially encouraging is the fact that the PTZR 
approach provides such consistent results without the need of collecting gauge data; a time 
consuming and costly effort. Gauge data are absolutely necessary for validation of results, but in 
a regime such as Kwajalein with limited locations for gauge sites, and logistically difficult 
maintenance and data collection, the PTZR approach can provide high quality near-real-time 
radar rain estimates when QC’d and calibrated radar data are available. 

5. Summary 

DP radars are a vital tool for GPM validation due to their applications to rainfall 
microphysical retrievals. The ability to provide consistent and long-term calibrated ground- 
based DP measurements will prove essential for calibration of the core GPM satellite and for 
development of physically based passive microwave radiometer algorithms over land. Through 
the retrieval of DSD parameters relating drop size and shape, rainfall estimation, and 
hydrometeor identification, the DP radar can provide validation of parameterized microphysical 
properties. 

The development and adaptation of algorithms for QC, calibration, rain rate estimation, 
and hydrometeor identification to be applied in a continuous, operational environment has been 
initiated with DP measurements from KPOL in an oceanic environment. Extension of these 
applications to ground-based WSR-88D radars is expected to begin in late 2010, dependent upon 
the polarimetric upgrade schedule. Presented are operational algorithms for QC, and absolute 
reflectivity calibration using polarimetric properties of the rain medium. Application of the QC 
algorithm has shown to be robust with superior results compared to the standard TRMM GV QC 
algorithm that employs height and reflectivity thresholds. The ability to detect and remove 
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ground and sea clutter embedded in precipitation echo is a distinct advantage of the DP 
algorithm. Through application of thresholding tests for correlation and standard deviation of 
differential phase, almost all clutter-type returns are identified and removed. In contrast, the 
TRMM GV algorithm can remove ground clutter only when not embedded in precipitation, and 
can be a labor-intensive process. In addition, the RVP8 processor correctly identifies multiple- 
trip echo through an automated SNR power test; a successful result to which the QC algorithm 
takes full advantage. 

A technique to determine Zdr calibration through analysis of combined Zdr Zh profiles 
was developed and applied to KPOL data. This application relies on the independent RCA for 
reflectivity calibration, and adjusts Zdr profiles to match a disdrometer distribution when 
birdbath scans are not available. By this technique, uncertainty has been mitigated in Zdr data 
from significant rainfall events in 2006 and 2007 and has allowed application of self-consistency 
reflectivity calibration. 

A published phase-based self-consistency approach to determine absolute reflectivity 
calibration using properties of the rain medium in light rain has been tested with KPOL data 
from five case studies and found to provide good results (within ± 1.0 dB) as compared to the 
independent RCA method. The approach follows the work of Ryzhkov et al. 2005 where K D p 
data from light rain events are integrated and compared against a model consistency equation 
derived from disdrometer data at Kwajalein. In lower reflectivity bins (30-35 dBZ), the observed 
Kdp measurements indicate a high bias relative to the disdrometer-based model, possibly due to 
difficulty in extracting the true Kdp signal. In mid-reflectivity bins (36-43 dBZ), there is good 
agreement between observed and model Kdp profiles (with one exception). The results indicate 
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that the method can be successfully applied to lighter precipitation regimes. The most consistent 
results were found in cases with uniform rain shields. 

The benefit of QC’d and calibrated radar data was demonstrated via rainfall estimation 
using the polarimetrically-tuned Z-R approach of Bringi et al. (2004). With corrected radar data, 
rainfall estimates showed very good agreement with independent rain gauges. The absolute bias 
between radar and gauge accumulations was reduced from 74% to near 16% when QC’d and 
calibrated radar data were used. This result is impressive considering that the Bringi approach 
does not require gauge data for calibration and is therefore able to provide high-quality radar rain 
products in near-real-time. Ongoing work will include quantitative rainfall estimation, QC aloft 
(melting zone and above) and full volume qualitative hydrometcor identification. 
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List of Figure Captions 


Fig 1. The location of Kwajalein Atoll in the Republic of the Marshall Islands (from Wolff et al. 
2005) The KPOL S-band radar is located on Kwajalein Atoll at the center of the image (8.7°N, 
167.7°E). Rain gauge locations from the KWA network are shown as black squares. 


Fig. 2. Comparison of co-polar pnv (top panel), Zdr deviation (Zdr - average Zdr) (middle), and 
Kdp (bottom) polarimetric parameters in light rain from CHILL, SPOL, and KPOL radars. 
CHILL observations are from 16 July 2004 2142 UTC, SPOL observations are from 26 Jan. 
1999 2139UTC, and KPOL observations are from 1 1 Aug. 2007 0352 UTC. 


Fig. 3. KPOL quality control algorithm flowchart for the first two elevation angles. All data are 
below the melting level. 


Fig. 4. QC examples showing persistent ground clutter (a) in raw KPOL data along atoll edges 
(50 km range ring), and removal of embedded and non-embedded ground clutter by the QC 
algorithm (b). A raw 0.4 degree elevation sweep (c) shows significant multiple-trip echo 
(southwest of KPOL, 20-160 km range), sea clutter, and typical noise. The QC’d field (d) shows 
dramatic improvement. Range rings are at 50 km intervals. 


Fig. 5. Map of the clutter field at Kwajalein. Range rings arc drawn at 10-km intervals from the 
radar site. Figure from Silberstein et al. 2008. 
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Fig. 6 Distribution of average Zdr measurements by reflectivity bin. The unadjusted Zdr curves 
(top panel) show the bias in five case studies relative to the JWD Zdr reference (bold, no-symbol 
line). The bottom panel shows Zdr distributions bias-adjusted to the JWD reference. Zh 
distributions have been calibration corrected via the independent RCA method. 

Fig. 7. Measurements of reflectivity with height (top left), mean Zdr with height (top right) 
correlation with height (bottom left), and the mean Zdr profile with azimuth (bottom right) for a 
shallow light rain event at Kwajalein from 07 Mar 2008. The correlation profile indicates 
sampling in the rain medium at the 3.0 km height. The azimuth profile shows the influence of 
directional-dependent response. The mean Zdr bias of +1.12 dB (heavy dashed lined) was 
determined over the full azimuth cycle. The dash-dot lines represent one standard deviation. 

Fig. 8. Comparison of radar processor determined Kdp measurements and computed Kdp 
measurements from self-consistency theory. The self-consistency relationship was developed 
from disdrometer measurements at Kwajalein. Results are shown from each of the five case 
studies. 
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Figure 9: Scatterplots of independent gauge and radar monthly accumulations for the period 
July-December 2008 at Kwajalein using the polarimetrically-tuned Z-R approach of Bringi et al. 
(2004). Results show the benefit of quality controlled and calibrated radar data (left panel) vs. 
uncorrected data (right panel). The absolute bias using corrected and uncorrected radar data are 


16.5% and 74.5%, respectively, with bias calculated as 


Radar -Gauge\ . 

. 1 . Regression lines 

Gauge 


(solid) are shown for each dataset along with the 1-to-l line. 
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Fig. 1 : The location of Kwajalein Atoll in the Republic of the Marshall Islands (from Wolff et al. 
2005) The KPOL S-band radar is located on Kwajalein Atoll at the center of the image (8.7°N, 
167.7°E). Rain gauge locations from the KWA network are shown as black squares. 


37 







Fig. 2: Comparison of co-polar p H v (top panel), Z DR deviation (Z DR - average Z DR ) (middle), and 
Kdp (bottom) polarimetric parameters in light rain from CHILL, SPOL, and KPOL radars. 
CHILL observations are from 16 July 2004 2142 UTC, SPOL observations are from 26 Jan. 

1999 2139UTC, and KPOL observations are from 1 1 Aug. 2007 0352 UTC. 
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Fig. 3: KPOL quality control algorithm flowchart for the first two elevation angles. All data are 
below the melting level. 
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KWAJ 05 Apr 2006 06:38:14 UTC CZ Elev: 0.41 



NASA/TRMM Office 

KWAJ 05 Apr 2006 04:28:08 UTC CZ Elev: 0.41 



Fig. 4: QC examples showing persistent ground clutter (a) in raw KPOL data along atoll edges 
(50 km range ring), and removal of embedded and non-embedded ground clutter by the QC 
algorithm (b). A raw 0.4 degree elevation sweep (c) shows significant multiple-trip echo 
(southwest of KPOL, 20-160 km range), ground and sea clutter, and typical noise. The QC’d 
field (d) shows dramatic improvement. Range rings are at 50 km intervals. 
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Fig. 5: Map of the clutter field at Kwajalein. Range rings arc drawn at 10-km intervals from the 
radar site. Figure from Silberstein et al. 2008. 
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Reflectivity (Zh) [dBZ] 

Fig. 6: Distribution of average Zdr measurements by reflectivity bin. The unadjusted Zdr curves 
(top panel) show the bias in five case studies relative to the JWD Zdr reference (bold, no-symbol 
line). The bottom panel shows Zdr distributions bias-adjusted to the JWD reference. Zh 
distributions have been calibration corrected via the independent RCA method. 
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Fig. 7: Measurements of reflectivity with height (top left), mean Zdr with height (top right) 
correlation with height (bottom left), and the mean Zdr profile with azimuth (bottom right) for a 
shallow light rain event at Kwajalein from 07 Mar 2008. The correlation profile indicates 
sampling in the rain medium at the 3.0 km height. The azimuth profile shows the influence of 
directional-dependent response. The mean Zdr bias of +1.12 dB (heavy dashed lined) was 
determined over the full azimuth cycle. The dash-dot lines represent one standard deviation. 
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Fig. 8: Comparison of radar processor determined Kdp measurements and computed Kdp 
measurements from self-consistency theory. The self-consistency relationship was developed 
from disdrometer measurements at Kwajalein. Results are shown from each of the five case 
studies. 
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Figure 9: Scatterplots of independent gauge and radar monthly accumulations for the period 
July-Decembcr 2008 at Kwajalein using the polarimetrically-tuncd Z-R approach of Bringi et al. 
(2004). Results show the benefit of quality controlled and calibrated radar data (left panel) vs. 
uncorrected data (right panel). The absolute bias using corrected and uncorrected radar data arc 


16.5% and 74.5%, respectively, with bias calculated as 


Radar - Gauge | 
Gauge 


Regression lines 


(solid) are shown for each dataset along with the 1 -to- 1 line. 
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Tabic 1. KPOL Basic Characteristics, Moments Recorded and Field Descriptions 


Processor: SIGMET RVP8 / RCP8 combination 

Frequency: 2.8 GHz 

Wavelength: 10.71 cm 

Beamwidth: 1.1° (horizontal and vertical) 

Operation Mode: Linear, horizontal and vertical 

Simultaneous dual transmit and receive (STAR) 

Moments Recorded: 

Zh: Reflectivity, horizontal component 

V r : Radial velocity 

o v : Spectrum width 

Zdr: Differential reflectivity 

Odp: Total differential phase 

Kdp: Specific differential phase (processor computed) 

Phv( 0): co-polar cross correlation at zero lag 

Field Variable Descriptions [with associated units] for KPOL data: 

ZT [dBZ]: Raw reflectivity, horizontal component 

DZ [dBZ]: SIGMET ground-clutter corrected reflectivity, horizontal component 

CZ*[dBZ]: Quality controlled reflectivity, horizontal component 

VR [m s' 1 ]: Radial velocity 

DR [dB]: Differential reflectivity 

PH [°]: Total differential phase 

KD [° km' 1 ]: Specific differential phase 

RH [no units]: Co-polar cross correlation 

SW [m s' 1 ]: Spectrum width 

SQ [no units]: Signal Quality Index 

* The CZ field is added by the NASA TRMM Satellite Validation Office and contains the 
reflectivity data that have been corrected for non-precipitation echo and absolute calibration 
error. 
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Tabic 2. Task configuration of the KPOL radar. Columns arc task name, radar polarization, 
elevation angles (deg), and pulse repetition frequency (PRF). Volume scans alternate between A 
and B, with one surveillance scan between volume scan sets. There are 10 volume scans per 
hour (5-A scans and 5-B scans). The alternating scanning was replaced with a single 17- 
elevation volume scan with 6 volume scans per hour in October 2008. 


Task Polarization 


Elevation Angles (deg) 


PRF 


GVVOLA Dual 
GVVOLB Dual 
Surv TRMM Horizontal 


0.4, 1.4, 2.3, 4.2, 6.1, 8.0, 9.9, 11.8, 14.0, 16.6, 19.6,23.2 960 
0.4, 1.4, 3.3, 5.2, 7.1, 9.0, 10.9, 12.9, 15.2, 18.0, 21.3, 25.3 960 
0.4 396 


Volume Scan Strategy beginning October 2008 (Alternating volumes A and B replaced) 

GVVOL Dual 0.4,1.5,2.6,3.7,4.8,5.9,7.0,8.1,9.2,10.3,11.6,13.2, 

15.2, 17.7, 20.8, 24.6, 29.2, 34.7 960 
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Table 3 Empirical investigation of relative DP radar performance in light rain. For both CHILL and SPOL radars, one case was 
selected for comparison. For KPOL, 12 cases from 2006 and 2007 were analyzed. The KPOL statistics represent the range of values 
observed from these cases. The standard deviation of the measured differential phase o(Odp) was computed at each range gate from a 
running, centered 21 gate sample for CHILL and SPOL, and a 25 gate sample for KPOL corresponding to the length scale used for 
Kdp calculation. The average absolute deviation of Zdr (AAD Zdr- see equation 2) represents an upper-end estimate of the random 
measurement error. 


Radar 

Name 

Date/Time (Z) 
YYMMDDHHMM 

# of Gates 

Median Zh 
dBZ 

Median Kdp 
deg km" 1 

Median pnv 

Median o(ODP) 
deg 

AAD Zdr 
dB 

CHILL 

0407162142 

4946 

22.6 

0.0 

0.99 

2.5 

0.4 

SPOL 

9901262139 

6827 

22.5 

0.0 

0.99 

1.4 

0.3 

KPOL 

1 2 cases 

3900-25000 

20-28 

0.0 

>0.98 

2.2-3. 1 

0. 3-0.4 



Tabic 4. Case studies and time periods selected (HHMM UTC) for application of the DP QC 
algorithm. Statistics represent the maximum percentage of gates flagged as “bad” or “noise” by 
the individual QC tests from the first two elevation angles within each volume scan. Statistics 
are from the rainfall estimation range of 0-100 km. Gates flagged in total reflectivity (ZT) 
includes those from the pnv, ^(^dp), low reflectivity threshold, and missing gate tests from the 
DP fields. The ZT percentages vary significantly from case-to-case and arc heavily influenced 
by the consideration of gates less than 5 dBZ as “noise”. 

Case 

pHV 

KDP 

ZDR 

o(ODP) 

Total ZT 

19 Dec 2006 
0544-0637 UTC 

3 

1 

17 

2 

5 

03 Aug 2007 
0200-0453 UTC 

2 

2 

34 

3 

8 

1 1 Aug 2007 
0325-0439 UTC 

8 

3 

30 

6 

29 

19 Sep 2007 
1542-1729 UTC 

4 

2 

33 

4 

24 

23 Nov 2007 
1603-1730 UTC 

5 

2 

23 

7 

29 



Tabic 5. Comparison of DP and non-polarimetric TRMM GVS QC results with respect to 
ground clutter detection. Raw reflectivity data were extracted from locations with frequent 
ground clutter echoes. The gate count represents the remaining gates considered to be ground 
clutter after QC is applied. For this study, gates > 55 dBZ are considered ground clutter. DP QC 
shows superior performance as compared to the non-polarimetric (TRMM GVS) QC as 
evidenced by the low gate counts. 


Case Number of unedited gates DP QC GVS QC 

> 55 dBZ / total extracted gates Gate Count Gate Count 

>55 dBZ >60 >65 >55 dBZ >60 >65 


19 Dec 2006 

71 / 11907 

0 

0 

0 

53 

3 

0 

03 Aug 2007 

237/35721 

3 

0 

0 

193 

45 

0 

1 1 Aug 2007 

383 / 

15720 

0 

0 

0 

85 

31 

0 

19 Sep 2007 

137/ 

18522 

0 

0 

0 

111 

28 

0 

23 Nov 2007 

528/ 

18429 

9 

1 

0 

478 

166 

34 
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Table 6. Reflectivity bin statistics to determine Kdp data validity from the 19 Dec 2006 case 
study. Shown are the reflectivity bin, average Kdp within the bin, number of independent Kdp 
samples, times reduction in Kdp standard error due to averaging, maximum allowed standard 
deviation of Kdp, and the standard deviation range of Kdp within the bin. 


dBZ Bin 

Kdp Avg 

(° km' 1 ) 

#indep Kdp 
Samples 

Times 

Reduction 

Max allowed 
o(Kdp) (° km' 1 ) 

Observed Data Range 
o(Kdp) (° km’ 1 ) 

30.0 

0.053 

1741 

41 

0.011 

0.004 - 0.006 

31.0 

0.059 

1833 

42 

0.012 

0.004 - 0.006 

32.0 

0.066 

1915 

43 

0.013 

0.004 - 0.005 

33.0 

0.076 

1957 

44 

0.015 

0.004 - 0.005 

34.0 

0.088 

1906 

43 

0.018 

0.004 - 0.005 

35.0 

0.102 

1790 

42 

0.020 

0.004 - 0.006 

36.0 

0.120 

1615 

40 

0.024 

0.004 - 0.006 

37.0 

0.141 

1402 

37 

0.028 

0.005 - 0.006 

38.0 

0.170 

1158 

34 

0.034 

0.005 - 0.007 

39.0 

0.210 

945 

30 

0.042 

0.006 - 0.008 

40.0 

0.260 

763 

27 

0.052 

0.006 - 0.009 

41.0 

0.321 

616 

24 

0.064 

0.007-0.010 

42.0 

0.391 

501 

22 

0.078 

0.008-0.011 

43.0 

0.457 

402 

20 

0.091 

0.008-0.012 

44.0 

0.543 

309 

17 

0.109 

0.010-0.014 

45.0 

0.626 

236 

15 

0.125 

0.011 -0.016 

46.0 

0.708 

176 

14 

0.142 

0.013-0.018 

47.0 

0.822 

127 

12 

0.164 

0.015-0.021 

48.0 

0.913 

79 

10 

0.183 

0.019-0.027 
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Table 7. Self-consistency derived calibration adjustments as compared with the independent 
RCA approach, and the absolute value of their difference. Shown is the cumulative number of 
Kdp samples compared (measured and derived via consistency) from the reflectivity bins in the 
Zmin (30 dBZ) to Zmax (48 dBZ) range. 


Case 

Kdp 

samples 

Self-Consistency Calibration 

ZeiAs(dB) 

RCA 

(dB) 

| Difference 
(dB) 

19 Dec 2006 

546434 

-2.44 

-1.95 

0.49 

03 Aug 2007 

1521223 

-2.06 

-1.78 

0.28 

1 1 Aug 2007 

335658 

-1.46 

-1.91 

0.45 

19 Sep 2007 

400012 

-0.93 

-1.91 

0.98 

23 Nov 2007 

740603 

-2.17 

-2.46 

0.29 
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